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Abstract 



Mass dependence of the effective source radii, observed in hadronic 
■ Z° decays by several LEP I experiments, is analyzed in a model which 

assumes proportionality between four-momentum of a produced particle 
Q_i| and the four- vector describing its space-time position at the freeze-out. 

O ■ It is shown that this relation (commonly accepted in description of high- 

energy collisions) can explain the data, provided all particles are emitted 
from a "tube" of ~ 1 fm in diameter at a constant proper time ~ 1.5 fm. 



1 Introduction 



Recently, two of us have pointed out [1] that the generalized Bj or ken-Gottfried 
hypothesis [2], relating the space-time position of a hadron produced in a high- 
energy collision to its 4-momentum, can qualitatively explain the mass-dependence 
of the interaction radii observed in e + e~ annihilation at LEP I [3, 4]. As discussed 
in detail in [1], this effect is a manifestation of the well-established observation 
[5] that a correlation between the momentum and the emission point of a parti- 
cle can drastically affect the results of the HBT correlation experiment. In the 
present note we want to explore in more detail the idea formulated in [1] in order 
to verify, if it indeed provides a viable framework for the understanding of the 
mass dependence of the HBT parameters. 
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The generalized Bjorken-Gottfried hypothesis, as formulated in [1], postulates 
the linear relation between the 4-momentum of the produced particle and the 
space-time position at which it is produced 1 : 

g M = Ax M . (1) 

Relation (1) implies 



where M\ = E 2 — q 2 = m 2 + qj_, and 

T = VW^z 2 (3) 

is the longitudinal proper time after the collision (t and z are time and longitudinal 
position of the production point). 

Since this picture is purely classical, it represents only a qualitative idea, 
whose application to the description of the actual data requires an adequate re- 
formulation taking into account the effects of the quantum nature of the system 
considered. In [1] we have proposed to use (1) and (2) as a guide-line for con- 
struction of the source function 2 S(P, X) [7, 8] related to the density matrix in 
momentum space by a Fourier transform 

p(q = P + \Q, q' = P - l -Q) = J d*Xe^ x S(P, X). (4) 

All variables are four dimensional, so that both space and time integrations are 
involved. Thus specifying the source function fixes completely the single particle 
properties of the model. 

Similarly as the standard Wigner function [9], S(P, X) gives approximately (as 
far as possible without contradicting quantum mechanics) the single-particle dis- 
tribution in momentum and in space-time. Therefore it has an intuitive meaning 3 , 
which can be exploited for implementation of the relations (l)-(2). 

The construction of S(P,X) (described in the next section) requires spec- 
ification of several details which can only be determined from a confrontation 
with data. We thus analyse the single- and two-pion distributions in the frame- 
work of the scheme proposed in [1] and compare them with the data from the 
DELPHI experiment at LEP I [10] in order (i) to pin down the parameters de- 
scribing the system and (ii) to verify that their values are reasonable, i.e., that 

1 To our knowledge, the first application of relation (1) to a discussion of the quantum 
interference between identical particles was proposed (in a different context) by Csorgo and 
Zimanyi [6] 

2 It was called there a "generalized Wigner function". 

3 It should be realized that, in contrast to the standard Wigner function which relates the 
particle wave functions at different positions but at the same time, the source function relates 
the particle production amplitudes at different positions and at different times, (as is clearly 
seen from (4)). Consequently some care is needed in order to assess its physical interpretation. 
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(2) 



the whole scheme is not unrealistic. Our conclusion is rather optimistic: the pro- 
posed scheme seems to describe adequately the LEP I data. As, in addition, it 
provides a reasonable intuitive picture of the hadronization, one may hope that 
it represents a useful framework for the description of the general properties of 
this process. 

In the next two sections we describe the construction of the source function 
and calculate the single particle density matrix in momentum space. In Section 4 
the formulae for two-particle correlation functions are derived. Comparison with 
the data is discussed in Section 5 and 6. The last section contains our conclusions 
and outlook. 

2 The Source Function 

To implement the conditions (l)-(2) we postulate the source function in the fac- 
torized form 

S(P,X)=F(t)S {1 S ± (5) 

where 



and 




Here 

X ± = t±z; P ± = P ±P Z . (8) 

so that 

Ml = P+P_; r 2 = X + X_. (9) 

We have used Gaussian forms in order to simplify the evaluation of the Fourier 
transform (4). As shown below, this admittedly crude assumption seems sufficient 
at the present level of analysis and accuracy of the data. 

The first exponential factor in S± represents a standard 4 cylindrically sym- 
metric "tube" of radius r± in configuration space 5 . The remaining exponential 
introduces a correlation between the momentum and the emission point of the 
particle, as required by the generalized Bjorken- Gottfried condition (1). Such 
correlations are known to influence strongly the HBT effect on particle spectra 

4 For related approaches see, e.g., [11]. 

5 To simplify the argument, we ignore the rapidity and z dependence of the single particle 
spectrum. This seems a reasonable approximation in the central rapidity region at high energy 
and can be easily removed, if necessary. 
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[5] and were shown in [1] to induce the mass dependence of the HBT radii. The 
reason for the mass dependence is easily seen from (2): at fixed r, the correlation 
between momentum and position depends on transverse mass. 

The parameters S\\ and 8± parametrize the correlation lengths or - in other 
words - the size of the region (centered at X M ) from which the observed particles 
emerge. They are arbitrary, apart from inequalities required for consistency with 
the quantum uncertainty principle, see, e.g., [14], which put on them some weak 
lower limits. 

Finally, the function F{r) gives the distribution of the proper time r at which 
the particles are created. 

Using now the formula (4) and the formulae of this section one can evaluate 
the single-particle density matrix in momentum space, which is a necessary step 
in the evaluation of both the single-particle spectrum and the two-particle HBT 
correlations. This calculation is outlined in the next section. 



3 Density matrix in momentum space 

Substituting (5) into (4) we have 

P(q,q') = J TdrF(r)p\ip ± 

where 



P± 



and 



P\\ = J dr ie V - 



Here 77 is the pseudorapidity: 



and 



Mi 



M 



X+) )+i(Q t-Q llZ ) 



T 



From (4) we also see that 



P = \(q + q'); Q = q-q' 



(10) 

(11) 
(12) 

(13) 
(14) 

(15) 



The Gaussian integral (11) for p± can be explicitly performed with the result 

.M±tv 2 - A 



P±(q±,q'±) = 2vrr e// exp 



P 2 Qjrjff 
2lu 2 2 



cxp 



-P±Qj 



(16) 
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where 

r 2 I? 

2 2 i s2 2 2 / 2 2 ' ± u ± /i r?\ 

u = M ± v +5 ± ; v = tJt ; r eff = — — (17) 

The longitudinal integral is somewhat more complicated. We first express the 
energies and longitudinal momenta in terms of the respective rapidities 

P± = M±e ±Y ; Q = m± cosh y — m! L cosh y'\ Q\\ — m± sinh y — m! L sinh y' 

(18) 

Substituting this into (14) we obtain after some algebra 

M? 

V = _A(i-cosh(y-77)) 

d ll 

(P±Q± cosh(F -r])+ m ± m' ± sinh(y - y') sinh(?? - F)) (19) 

One sees from this formula that V depends only on Y—rj and y—y'- Consequently, 
after integration over r\ the result depends only on y — y' . 

The integral (12) can now be evaluated. We first change the variable 77 — F — > 
£. Then we express the hyperbolic functions of £ by exponentials. Using the 
integral representation 

POO 

K {z) = / e~ zcoshy dy (20) 
Jo 



we obtain 



M 2 

p || =2expl 7 f ItfoOO (21) 



with 



A/f 4 t 2 — » — » rMi - -* t 2 

s- = - ^r(P ± Qx) 2 - 22— ^P ± g ± + —jmlm'lsmh^y - y') (22) 



<5f M 2 ^ 5\ ^ AM 2 

Using the identities 



cosh y + m' ± cosh y' = 2Mj_ cosh F 

sinh y + m' ± sinh y' = 2M± sinh F (23) 
we arrive, after some algebra, at our final formula 

s 2 = ^-r 2 Q 2 t - i T -^(ml-m'l) (24) 

where 

Qt =Ql-Q\ (25) 

This completes the evaluation of the single-particle density matrix. 
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4 Single-particle distribution and two-particle 
correlation function 



The single particle distribution is given by the diagonal elements of the density 
matrix. From the formulae of the previous section we thus obtain 

^^ = 2 ^ exp (f)^(f) /tei) (26) 

where 

I(ql) = J TdrF(r)Q- 2 exp (-|L)) (27) 

and 

to 2 = m 2 ± v 2 + S 2 ± (28) 

To obtain information on the two-particle correlation function, one has to 
make further assumptions. We follow the standard treatment [7, 8, 14], assuming 
that one can evaluate the two-particle correlation function as if there were no 
other correlations between particles except for those induced by quantum inter- 
ference. Under this condition the normalized two-particle correlation function is 
given by 

The final point we want to discuss is the selection of the variables used for com- 
parison with the data, and the corresponding Jacobians. 

The system of two identical particles is fully described by the 6 components 
of the momenta. The phase-space volume is 

dtt = dydy f d 2 q ± d 2 q' ± = 2ndydy'dPldQ 2 ± d(f) (30) 

— * — * 

where <p is the angle between P± and Q±. 

As seen from the formulae of Section 3, it is convenient to replace the two 
rapidities by Q\ and Y. Using the identity 

Qt = ( m i-L - r m 2 ±) 2 - 4m 1± m 2 _L sinh 2 ^ — - V2 ^j (31) 



one obtains 



d Vl dy 2 = . d lM (32) 

yjQl - 2(m 2 ± + m 2 ± )Q 2 + (m 2 ± - m 2 ~* 



L 2±J 



Thus everything can be expressed in terms of Y, Q 2 , Q 2 ± , P 2 and 0. It is sometimes 
convenient, however, to consider other variables. For completeness we give below 
some kinematic relations which can be useful in data analysis. 

/i 2 = 2(m 2 1± + m 2 2± ) = 4m 2 + 4P 2 + Q\ (33) 
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Qout = \Q±\ cos0; Mi = /i 2 -g t 2 
m l± ~ m i±_ = 2\P±\\Q±\ cos0 = 2\P±\Q 0Ut 



(34) 
(35) 



5 Data: Single-particle distribution 

Eq.(lO) represents the density matrix as an integral over the proper time r at 
which the particles are produced. In the present paper, following [1], we shall 
accept the approximation that the production happens in a very narrow interval 
of r, so that the integration over r simply amounts to introduce a fixed value 
r = r in the formulae of the previous section. In this way the unknown function 
F(t) is replaced by one parameter, r . The other three parameters are: v, <5j_ 
and <5||, each with a very clear physical meaning. 

The first step in the data analysis should be the correct description of the 
single particle spectrum using the Eq.(26). This also allows to restrict somewhat 
the values of the four parameters we have to our disposal. 

For this purpose the data sample of ~ 3 * 10 5 Z° hadronic decays from the 
DELPHI experiment was used [10]. Hadronic events have been selected using 
standard DELPHI cuts [12] which limit the contamination of beam-gas, 77,/x/i 
and ee events to less than 0.1% and rr contamination to below 0.2%. In this 
study only two-jet like events were considered (Thrust > 0.95). The distribution 
of transverse momentum with respect to the event thrust axis was constructed 
for tracks which are not obvious decay products of V 0, s and are not positively 
identified as K's or p's. 

The formula (26) has been fitted to the distribution in the q\ region up to 
2.0 GeV 2 which contains 99% of the data. The small q\ < 0.020GeV 2 interval 
was excluded from fits because it is particularly uncertain experimentally due to 
the contribution of low momentum pions from the D* decays and to a possible 
distortion of the distribution by imperfect approximation of the primary quark 
direction by the measured thrust axis. 

The unknown parameters (v, 5± and 5u plus normalization) were determined 
by minimizing a x 2 function with the MINUIT program [13]. The parameter tq 
was kept fixed in fits at the value of 1.5 fm. Since correlations between some fit 
parameters are high, we used a 3 standard deviation correlation ellipse for each 
parameter pair to estimate their statistical errors. In addition, the stability of 
the results of the fit has been checked by varying the range of fitted q\. The 
extreme values of a parameter resulting from these fits were used as a measure 
of the range of a parameter allowed by the data. 

Although the fits converged easily to a deep minimum, their x 2 value was 
not always statistically acceptable. Nevertheless, we accept these fits as long as 
the shape of the experimental distribution is correctly reproduced over a wide 
range of q\ . The too high Xmin/dof value is attributed to underestimated errors 
(imperfect treatment of systematics, no background subtraction etc) rather, than 
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to failure of the fitted hypothesis. An example fit is shown in Fig.l. The fit results 
are summarized in Table 1. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

q ± 2 [GeV 2 ] 



Figure 1: The fit of the formula (26) to q\ distribution of pions, (v = 0.286 , 
S ± = 0.303 , 5|| = 0.172) 



6 Two-particle correlation function 

The two-particle correlation function is calculated according to formula (29) us- 
ing the values of the parameters from Table 1. Correlations between parameters 
determined by the fit are taken into account by sampling the area of the 3a 
correlation ellipse for each parameter pair with a Gaussian distribution. At the 
presence of correlations between the parameters this procedure results in a much 
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Table 1: Model parameters determined in the fit. 



V 


S ± { GeV ) 


S\\{ GeV ) 


X 2 /dof 


o.286±8:835 


n ^n3 +a010 




397/193 



more reliable estimate of the statistical error on the calculated correlation func- 
tion value than the 1-dim. errors given by MINUIT. Systematic effects of the fit 
are accounted for by repeating calculations for extreme values of the parameters 
allowed by the q\ distribution. This procedure was needed especially for the 
parameter 5jj which is the least constrained by the one-particle distribution. In 
this way the value of the correlation function for each particle was determined 
together with its error, where both error contributions described above have been 
added linearly. The result of the calculation for pions is shown in Figs.2a,b for 
the perpendicular and parallel components, respectively. 




0.2 0.4 0.6 0.8 1 1.2 1.4 0.2 0.4 0.6 0.8 1 1.2 1.4 

Q ± 2 [fin" 2 ] Q||2 [fin 2 ] 

(a) (b) 



Figure 2: Correlation function vs Q\ (a) and Q?, (b) for it 
One sees that these functions are not gaussian. Nevertheless, to compare 
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with the existing data, we have approximated the results by Gaussians in the 
region O.lOOGeV < Q < 0.250Ge\/ (Fig. 2). The exclusion of the small Q 2 region 
may be justified by the well-know fact that the measurements of the correlation 
function at very small Q 2 are very uncertain (and often this region is omitted in 
the data analyses at LEP). The results of this procedure for pions are shown in 
Fig. 3(a) together with measurements of the three LEP experiments [3] done in 
correlation studies in two- and three dimensions. 




0.12 0.14 0.16 0.2 0.4 0.6 0.8 1 1.2 



m[GeV] m[GeV] 

(a) (b) 

Figure 3: R\\ and R± calculated from the model (shaded bands) : (a) for ir, (b) 
for 7r,K,p,A. Data points in (a) represent results of 2- and 3-dimensional analyses 
of LEP data [3]. Data points in (b) represent 1-dim source radius R [4] 

In the experimental studies Q 2 is decomposed into three components in the 
longitudinal center-of-mass (LCMS) frame, where the sum of the three-vector 
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momenta is perpendicular to the thrust axis. Choosing Ql parallel to the thrust 
axis, Q out parallel to the sum of the momenta (in the LCMS frame) and Q S ide 
orthogonal to both, the decomposition reads: 

1 2 

Q 2 = Qi+Qi de +Qi t a-p 2 y, P = Po ; t l P ; ut m 

Hi l + Cj2 

DELPHI performed the analysis in two dimensions, determining correlation radii 
corresponding to Q% = Q 2 side + Q 2 out ( R ± ) and to Q\ [R\\). The L3 and OPAL 
analyses where done in three components, therefore, as the measure of R± we have 

chosen R± = y — 2 Slde . The spread of the experimental results is considerable, 

which might be due to different methodology (eg. different reference sample, 
corrections for Coulomb repulsion etc) and differing phase space regions selected. 
However the clustering of the data points for R\\ at larger values than those for R± 
is clearly visible. Shaded bands in the plot represent the errors on the calculated 
R± and R\\. 

The most remarkable feature seen in this figure is that the model predicts, in 
agreement with experiment, a longitudinal radius much larger than the transverse 
one. One also sees that both radii fall in the ballpark of about 1 fm which is 
hardly surprising. It is also clear that -within the rather large theoretical and 
experimental errors- the model seems to account for the gross features of the 
data. 

The calculations were also done for kaons, protons and A's. Since the available 
data samples for these particles are much smaller than those for pions, their 
q\ distributions are not discriminative enough to pin-point reliably the model 
parameters. Therefore the parameters were taken the same as for pions (given 
in Table 1). This assumption is to be verified once better data are available but 
we have checked that it reproduces reasonably well the main characteristics of 
the transverse momentum distributions of kaons, protons and A's. The resulting 
correlation functions for all these particles are of Gaussian form over a wide range 
of Q\ and Q 2 (up to several hundred MeV) and thus the source radii are well 
defined. 

The mass dependence of the calculated R\\ and R± is plotted in in Fig. 3(b) 
where also the available results of the LEP experiments [4] are shown. 

One sees that the inequality R\\ > R± is still satisfied in the model, although 
the difference between the two radii at higher masses is not as large as in the 
case of pions. The data points in this figure correspond to the correlation ra- 
dius Rq determined in 1-dimensional analyses (the only available data for heavy 
particles) [4]. Multiple entries of the result from the same experiment for pi- 
ons correspond to different measurements made with different reference samples. 
The points at kaon mass represent measurements for both and K ± pairs. 
The measurements for A pairs come from spin analysis (except for the second 
ALEPH point with small error) where there is no need for a reference sample. 
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The correspondence between R and the two radii R\\ and R± is not obvious (at 
least experimentally), but the trend of the data is reasonably well reproduced 
by the model. More accurate data on kaons would be of great help to further 
elucidate this point. 

7 Conclusions and Comments 

In conclusion, we have found that the correlation between the momentum and 
the production point of a produced hadron, suggested by the Gottfried-Bjorken 
hypothesis of in-out cascade, seems to account (at least approximately) for the 
observed correlation between identical particles observed in e + e~ annihilation. 
Together with data on the single particle transverse momentum distribution, it 
predicts strong anisotropy of the two-pion correlation function, in agreement with 
the observations. The mass dependence of the "effective source radius" is also 
adequately described. Large uncertainties, both in the theoretical determination 
of the model parameters, and in the experimental data do not allow, however, to 
obtain more quantitative conclusions. 
Several comments are in order. 

(i) One sees from the Table 1 that out of the three parameters determined 
from transverse momentum distribution, only two (v and S±) are relatively well 
constrained, whereas 5\\ is only weakly restricted. It is interesting to note that 
- within the wide range allowed - the value of 5\\ is consistent with the relation 
5\\ = 5± corresponding to the quasi-isotropic case. Accepting this isotropy, we see 
that the effective width of the momentum distribution around the average given 
by the Bjorken- Gottfried condition (1) is close to 300 MeV, a rather reasonable 
value. 

The small value of v is also recomforting, as it guarantees that the transverse 
expansion of the "tube" from which the particles are emitted is not unreasonably 
fast. 

(ii) The mass dependence of the effective HBT correlation radii calculated 
from the model and shown in Fig. 6 was obtained under the assumption that the 
parameters of the model do not depend on particle masses. This assumption was 
verified (within large experimental uncertainties) for v, 5± and 5\\. No such direct 
check is available for r . However, since the observed mass dependence of the 
HBT correlation radii does agree -at least approximately - with the experimental 
observations, we may take it as an argument that also r does not depend on 
particle masses: 

r ~ const(M ± ). (37) 

This seems a rather non-trivial conclusion, as it indicates that -within the Bjorken- 
Gottfried hypothesis (1)- all particles are emitted at, roughly, the same proper 
time t . 
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This result may be contrasted with a simple expectation (sometimes justified 
by uncertainty principle) suggesting 



1 



(38) 



r ~ 



M 



which would give a stronger drop of R with increasing mass, particularly for 
transverse radius. 

It seems that also the models based on the picture of string decay [15] would 
not fulfill the condition (37) but rather obey 



Such a relation corresponds, in our model, to a much weaker dependence of the 
correlation radii on the particle masses and thus an other mechanism would have 
to be invented to describe the data [16]. 

(iii) In the present paper, studying the two-particle distribution, we considered 
-following the approach employed in experimental analyses [3]- the boost invariant 
variable Q\ and the variable Qf evaluated in the longitudinal center-of-mass 
system. Assuming boost invariance and azimuthal symmetry of the distributions, 
one finds that a complete analysis would involve 4 variables. As it is convenient 
to choose them boost invariant, one could use for instance the two transverse 
momenta |pu_| and |p 2 j_| the relative azimuthal angle 4>i ~ 02 and the relative 
rapidity y\ — y 2 . It would be interesting to see the data analysed in this way. 

(iv) Another interpretation of the experimentally observed HBT parameters 
was given in [17]. The authors take the point of view that the observed HBT radii 
do indeed correspond to the actual size of the particle emission region which is 
thus strongly dependent on the particle mass. They argue that this dependence 
may be understood from the uncertainty principle. It may be worth to point out 
that this approach is rather different from ours. In our description the parameters 
characterizing the particle emission region are mass independent and the observed 
change in the HBT radii comes solely from the momentum-position correlation 
as expressed in the assumed Bjorken-Gottfried condition (1). 

(v) Our prescription for construction of the two-particle correlation function 
out of the single-particle source function is similar to -but not identical with- that 
advocated recently by Heinz and collaborators [7, 18]. The difference is in the 
treatment of the fourth component of the average momentum vector of the two 
particles [P = (<fi + ^)/2]. As seen from (15), we use 



which guarantees the four-vector character of P M but brings it off mass-shell. 
Heinz et al. propose to take 



t ~ Mj_. 



(39) 



(40) 




(41) 
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which insures P M P M = m 2 but induces complicated transformation properties. 
In practice, however, one is only interested in the region q\ ps q 2 where these 
two prescriptions are not substantially different (the source functions themselves, 
however, are not the same). 
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